Fermi-Bose Mixtures Near Broad Interspecies Feshbach Resonances 
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In this Letter we have studied dressed bound states in Fermi-Bose mixtures near broad interspecies 
resonance, and implications on many-body correlations. We present the evidence for a first order 
phase transition between a mixture of Fermi gas and condensate, and a fully paired mixture where 
extended fermionic molecules occupy a single pairing channel instead of forming a molecular Fermi 
surface. We have further investigated the effect of Fermi surface dynamics, pair fluctuations and 
discussed the validity of our results. 



Since the observations of molecules of Fermi atoms 
near Feshbach resonances, fascinating pairing correla- 
tions in cold Fermi gases have been successfully inves- 
tigated both experimentally [lJ-Q and theoretically [sHU . 
Near broad resonances where the atom-molecule coupling 
is very strong, pair correlations can also be closely related 
to the ones in the BEC-BCS crossover theory pioneered 

a while ago num. 

Meanwhile, interspecies Feshbach res- 
onances in Fermi-Bose mixtures of 6 Li- 23 Na, 40 K- 87 Rb 
and 6 Li- s7 Rb have been experimentally observed 12- 16$. 
Recently, weakly bound 40 K- &7 Rb pairs prepared near 
Feshbach resonance were further successfully converted 
into cold molecules at JILA 17], which can potentially 
lead to exciting op por tunities of studying new quan- 
tum states of matter [18j. Previous theoretical studies on 
Fermi-Bose mixtures on the other hand have been mainly 
focused on narrow resonances or when the atom-molecule 
coupling is very weak |19l - l2l| ; phase boundaries in this 
limit depend on atom-molecule coupling strengths. 

Experimentally, creating and probing correlations in 
Fermi-Bose mixtures near narrow resonances are more 
challenging than near broad resonances. For a Feshbach 
resonance with a width AB and background scattering 
length flBG, an effective resonance energy width can be 



introduced as T res — 2m R a 2 BG (AfiAB) 2 /% [l(J,|22|; here 
Afi is the difference in magnet moments between the scat- 
tering and molecule channel and m R is the reduced mass 
of a pair of Bose and Fermi atoms. For interspecies res- 
onances so far confirmed [l2- 16, 23 1, r res can be a few 
orders of magnitude bigger than th e q uantum degener- 
acy energy at a density of 10 14 cto _3 [L2J. So some of well 
studied resonances are quite broad and many-body cor- 
relations near these interspecies resonances remain to be 
thoroughly understood. One of very fundamental ques- 
tions we hope to answer in this Letter is what is the 
nature of quantum matter near broad interspecies res- 
onances. Especially, how do dressed two-body bound 
states evolve when approaching resonances? and accord- 
ingly what kind of many-body correlations are developed 
in a quantum Fermi-Bose mixture? And physics near 
broad resonances can distinctly differ from that near nar- 
row resonances; some basic concepts introduced for nar- 
row resonances such as molecular Fermi surfaces might 
not be directly applicable here. Motivated by these con- 



siderations, we carry out a study on Fermi-Bose mixtures 
near broad interspecies Feshbach resonances which serves 
as a potential reference for more sophisticated analyses. 
Apart from the phase diagrams in terms of interspecies 
scattering lengths at,/ and Bose(ms)-Fermi(mf ) mass 
ratio ms/mF, we focus on bound state properties of 
near-resonance Fermi-Bose mixtures which can be po- 
tentially probed in experiments. Our results are useful 
for the understanding of correlations in 6 Li- 23 Na, 6 Li- 
s7 Rb, and 40 K- a7 Rb mixtures. For simplicity, we employ 
a simplest one-channel Hamiltonian which captures most 
important aspects near broad resonances, 
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where / k ,& k (/k, &k) are creation (annihilation) opera- 
tors for Fermi and Bose atoms respectively, and ^ = 
» fi ^ are kinetic energies for fermions (bosons) and f2 
is the volume. Vtf is the strength of interaction that is 
related to interspecies scattering lengths a^f via 
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here m R = tcibtcif / '{tub + itif), e k = ^ k 2 /2mfl. We 
assume that the background boson-boson interactions are 
are repulsive so that the mixture is stable; to illustrate 
the idea, here we only include interspecies scattering. 

Bound states and Pauli blocking effects We first 
consider the binding energy of a pair of Fermi and Bose 
atoms with opposite momenta (k, — k), in the presence 
of a condensate(BEC) and a Fermi surface of Fermi 
atoms which blocks all states below its Fermi momen- 
tum hkp- Pauli blocking effects of a Fermi sea indeed 
lead to dressed bound states at arbitrary small nega- 
tive scattering lengths but with an anomalous disper- 
sion, or a negative effective mass (see also discussions 
before EqH]). Furthermore, the energy Wb it takes to 
create a bound state from a non-interacting ground state 
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FIG. 1. (color online) The energy dispersion of bound states 
for mass ratio ms/jriF = 2.175 or 40 K- 87 Rb mixtures. From 
the top to bottom are Wb(Q,) in units of tp for 1/kpatf — 
-0.2,-0.1145,-0.05,0.0145. The shaded region represents 
pair excitation continuum. 



014 
0.3 






K-Rb — 
1 .i-Na 




B 0.2 


Li-Rb 


- 




\ \ 







i i i i ' 




-1.6 -1.2 -0.8 -0.4 
^^^__J/a bf k F 






1 3 6 9 12 15 

m B /m F 

FIG. 2. (color online) Scattering length a' 1 ' at which the 
energy cost Wb to create a molecule with HQ = from a 
Fermi gas-BEC mixture becomes zero, as a function of Bose- 
Fermi mass ratio ms/mF. Inset is Wb, the energy gap of 
molecules in units of ep as a function of kpa^f (and 1/at,/ < 
1/a' 1 ') for different mass ratios or mixtures. 



can be either positive or negative depending on scatter- 
ing lengths at/, a unique feature of Fermi-Bose systems. 
This is because to form a pair of atoms with opposite mo- 
menta (k, — k) near the Fermi surface |k| = kp, a Bose 
atom has to be promoted to right above the Fermi sur- 
face which results in an energy penalty ep — h 2 k F /2mB- 
For a bound state with an arbitrary total momentum Q 
or a kinetic energy £q = 7i 2 Q 2 /2(mp +mj), the energy 
cost is Wb{Q) = €p+€q+ujb- The Q-dependent binding 
energy u>b{< 0) can be obtained by solving the following 
two-body equation, 
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In the limit of small kpabf(< 0) and when Q = 0, 
Eq.© leads to uj b = -Aefexpij^j), e£ = h 2 k 2 F /2m R . 
The dispersion of bound states that can be probed us- 
ing photoassociative spectroscopy is shown in FigfT] [24| ■ 
For small negative scattering lengths a;,/, bound states 
are fully gapped with positive energies Wb(Q), and the 
ground state is a mixture of Fermi gas and BEC. How- 
ever, Wb(0), the energy gap of bound states vanishes 
at a critical scattering length a^ 1 ). In Figj2] we present 
results of versus ms/nip. For heavy Bose atoms, 
k F a^ approaches a small value of n/(\n\mp/4mB} +2). 

To ensure the stability of Q = molecules near the 
transition line we further examine M e ff, the effec- 
tive mass near Q = 0. At scattering lengths or when 
W B {0) = 0, we find 
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and the dimensionless function g(x 2 ) — x/[{l — 
„2\2n„ 1+x _i_ -2^-r)]. As far as m B /mp > 0.7 and the 
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energy penalty ep is not too heavy, M e // is positive, al- 
though it can be are much bigger than the bare mass 
tot(= tij? + tob) as a result of dressing in the Fermi sea. 
Below we focus on the limit of positive M e f f that is most 
relevant to the experimental mass ratio ms/mp (between 
2.175 and 14.5)[2J|. Note that the binding energy ujb is 
independent of the Bose atom density when the Fermi 
sea is treated as a static background. 

The above analysis at first sight seems to suggest that 
when Wb becomes negative, a small fraction of Fermi and 
Bose atoms start forming molecules or a dilute molecu- 
lar Fermi gas signifying a phase transition at a^. Such 
a picture was in fact previously proposed for mixtures 
near narrow resonances [l^, lipf . However since the ex- 
tent of molecules d m is typically comparable to or much 
longer than the Fermi wave length 2ix jkp near broad res- 
onances, pairs may be accommodated, even before the 
two-body gap Wb vanishes, in other more exotic forms 
without forming a molecular Fermi surface. Below we 
provide evidence for such a possibility. Note that a finite 
two-body gap Wb suggests a local stability of the Fermi 
gas-BEC mixture against emergence of a Fermi gas of 
molecules, i.e. when 1/a < \/a^\ a molecular Fermi gas 
can not be a ground state. It is in this limit that we illus- 
trate, based on an energetic analysis, that a third state 
or a fully paired state actually further lowers the energy. 

Energy landscape of pairing states Suggested by 
above discussions, we consider energetics of pairing states 
of Fermi-Bose atoms and from now on focus on homoge- 
neous mixtures with equal populations of fermions and 
bosons, i.e. Np — Nb- Although in principle pairing 
with a finite total momentum HQ can occur in ground 
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states, detailed calculations show that for a range of 
mass ratios (ms/nip > 0.2) that is relevant to Fermi- 
Bose mixtures so far studied in experiments, pairing in 
?iQ = channel is always dominating and favored, qual- 
itatively consistent with our analyses on the two-body 
bound states. Below we only show results of Q = pair- 
ing and adopt the simplest pairing wavefunction 



\g.s. >= exp(c bl) ]^(uk + Wk^-k + ^k/k)l uac > (5) 

fe#o 

where itk, Vk and ?7k are three families of variational pa- 
rameters. We obtain the energy of the variational states 
and then minimize it with respect to itk , ^k and r/k that 
are subject to the normalization condition |i<k| 2 + |fk| 2 + 
1 | 2 = 1. Equilibrium conditions can then be obtained 
and there are two solutions for any given k: i) a unpaired 
state with r;t = 1 and itk = t>k = 0; ii) a paired state 
with Vk = and vi = i (l - g), ul = i (i + g) , 

A = -^f Ek(l ~ »£)ukVk, where E k = y / (^+4A^ 
£k = e k ~ A 4 an d A 1 is the pair chemical potential. Pairing 
gap A, [i as well as the condensed population | Co | 2 are 

determined self-consistcntly, Np = Ek ( u k(l~ 7 7k)+ T /k) > 
A r B = |co| 2 + Ek«k(l-'7k)and 
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FIG. 3. (Color online), (a) a cr at which a first order phase 
transition occurs, as a function of mass ratio mj/niF. The 
dashed line is a' 1 - 1 along which a Fermi gas-BEC mixture be- 
comes locally unstable. Near a cr , there is a tiny window for 
phase separation (shaded area in the inset). Line 2 (1) is the 
shifted a cr (o^) due to quantum fluctuations (see discussions 
around Eq[7|. (b)E p , energy per pair of Fermi-Bose atoms in 
units of ep as a function of X, the size of Fermi surface of 
unpaired Fermi atoms for mB/mp — 2.175 or 40 K- 87 Rb mix- 
ture. From the top to bottom are E p for l/kpdbf = —5, 
-0.1145(for o cr ), -0.05, 0.0145(for a (1) ) and 0.2. X = 1 cor- 
responds to a unpaired state and X = is for a fully paired 
mixture. 
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To minimize the energy, we further choose rj^ to be a step 
function, -qy = 1 if |k| < Xkp, and zero otherwise; the 
dimensionless variational parameter X € [0, 1] specifies 
the size of residue Fermi surface of unpaired Fermi atoms. 
When Np = Ng, one can also verify that X 3 is equal to 
\cq\ 2 /Nb, i.e. the condensation fraction. 

In Fig |3l41 we present the main results of variational 
calculations. In FigOl the energy per pair of atoms 
is shown as a function of X, the size of Fermi surface 
of unpaired fermions. At small and negative scattering 
lengths, a Fermi gas-BEC mixture is a ground state and 
X = 1, fi = e| + W B (e£ = h 2 k F /2m F ) and A = 0. 
A fully paired state with X — becomes degenerate 
with the Fermi gas-BEC mixture (Xq = 1) at a criti- 
cal scattering length a cr , beyond which the paired state 
becomes a ground state. However, a Fermi gas-BEC mix- 
ture remains to be locally stable until scattering lengths 
reach the value of which is fully consistent with the 
above study of bound states. Our variational calcula- 
tions suggest a first order phase transition between a 
Fermi gas-BEC mixture and a fully paired mixture. This 
later state of extended molecules is conceptually differ- 
ent from a Fermi gas of molecules; instead, all molecules, 
though fcrmionic in nature, occupy the same pairing 
channel with zero total momentum. A direct compari- 
son of energies indicates that a fully paired mixture has 
lower energies than a Fermi gas of molecules provided 
l/a bf < l/a^K For i0 K- S7 Rb mixtures, l/k F a^ is 
about 0.25-0.3 and l/a^ > 1/aW; further towards the 
molecular side, a paired mixture is expected to evolve 
into a Fermi gas of molecules. 

In FigJU we further present the results on the pair 
breaking energy A and pair chemical potential [i. The 
pair breaking energy can be probed when applying rf 
pulses to transfer Fermi atoms to a different hyperfine 
spin state 0, [25| that weakly interacts with the Fermi- 
Bose mixture. The frequency shift in the rf spectroscopy 

should be ^Aw(k) = \ (f£ + +4A 2 ) . In the fully 

paired phase, Bose atoms are completely depleted and 
the Bose atom distribution ne(k) follows closely the 
Fermi atom distribution np(k). 

Fermi surface dynamics and pair fluctuations 
At small negative scattering lengths when quantum fluc- 
tuations are weak, our mean-field analyses on two-body 
bound states as well as the fully paired states are asymp- 
totically adequate. And a^\a cr estimated in the limit 
of large mass ratio (kpa^ ~ n/ (\n[mp / 4itib] +2) in 
Fig|3]) are quantitatively valid. However when kpa is of 
order of one, the fluctuations becomes substantial and 
the bound states can be further dressed in fluctuating 
particle-hole pairs; the energetic analyses are subject to 
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FIG. 4. (color online). (a)Pair-breaking energy A and (b)pair 
chemical potential /i in units of tp as a function of scatter- 
ing length kFdtf (near a cr ) for mass ratio mj/mp = 2.175 
or 40 K- 87 Rb. Dashed lines are for the values of meta-stable 
states. The tiny window for phase separation near a cr is too 
small to show here (see inset of Fig[3Ji). 



and and further stabilize the first order phase 

transition. Quantum Monte Carlo simulations similar 
to those in Ref. [29| remain to be carried out. 

In conclusion, we have examined dressed bound states 
and provided evidence of a new quantum state of ex- 
tended molecules near broad interspecies resonances. 
As far as three- and higher-body correlations are in- 
significant, our results can be applied to understand 
Fermi-Bose mixtures near broad resonances. We thank 
Gordon Baym, Immannuel Bloch, Kirk Madison, Subir 
Sachdev, Joseph Thywissen, Jun Ye and Zhenhua Yu 
for stimulating discussions. This work is supported by 
NSERC (Canada) and Canadian Institute for Advanced 
Research. Note added: Upon the submission of this 
work, we learned that finite-temperature mixtures near 
broad resonance were also studied in Ref. 30] . 



corrections. To clarify this, we estimate the dominating 
effects of Fermi surface dynamics, i.e. the Gorkov cor- 
rections ( GMB) in the two-body scattering vertex due to 
fluctuating particle-hole pairs [53, I2II 27], and the self- 



energy (SE) effect which mainly represents the effect of 
atomic Fermi surface smearing and the mass renormal- 
ization due to scattering by the condensate or Fermi sea. 
We then examine the pole structure of the propagator 
for a pair of Fermi-Bose atoms taking into account the 
vertex corrections and self-energy effects and obtain the 
binding energy ojb- The relative shift caused by the GMB 
effect (the leading order term R±) and the SE effect (the 
higher order term R2) is given as 



<Hn \u B \ 
In \ojb\ 



■ Ri(k F a b f) + i?2(fcFOb/) 2 In \k F a bf \... (7) 



where i?i ; 2 both are dimensionless quantities depend- 
ing on the mass ratio ms/rap- and can be obtained 
numerically [28]. For large mass ratios(ms/mi?), one 
finds that R\ ~ — \a.(m,B / Trip) / n an d Ri ~ — 2/(37r 2 ); 
and wb is reduced by a factor of mF/nriB solely due to the 
GMB effect. The net reduction in the binding energy ujb 
leads to a upward shift of l/a^ (line 1) in Fig|3] In addi- 
tion, we have estimated that, for the paired state, the am- 
plitude of pair fluctuations A\ ~ (A/ep) 4 ^ K/mjikp (n 
is the compressibility of the paired mixture), and the cor- 
responding zero point energy (in units of Cp) per particle 
A2 ~ (A/e^) 4 -y/m^fci?/K. We incorporate these quan- 



tum fluctuations (analogous to NSR effects [llj) which 
further favor pairing, and the GMB correction into our 
analysis of energetics of the paired state. For 40 K- S7 Rb 
mixtures, we find that A\ w 0.01, A2 « 0.003 and 
i?i w -0.4 and R 2 -0.1; the GMB and SE effects ap- 
pear to be more dominating. The modified critical lines 
obtained by extrapolating the above analyses to near res- 
onance are shown in Fig(3]and are qualitatively consistent 
with the mean field ones. Data suggest that quantum 
fluctuations tend to enlarge the window between l/a cr 
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